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ABSTRACT 

We report in this paper phase-resolved spectropolarimetric observations of the rapidly- 
rotating fully-convective M4 dwarf V374 Peg, on which a strong, mainly axisymmetric, 
large-scale poloidal magnetic field was recently detected. In addition to the original 
data set secured in 2005 August, we present here new data collected in 2005 September 
and 2006 August. 

From the rotational modulation of unpolarised line profiles, we conclude that 
star spots are present at the surface of the star, but their contrast and fractional 
coverage are much lower than those of non-fully convective active stars with similar 
rotation rate. Applying tomographic imaging on each set of circularly polarised profiles 
separately, we find that the large-scale magnetic topology is remarkably stable on a 
timescale of 1 yr; repeating the analysis on the complete data set suggests that the 
magnetic configuration is sheared by very weak differential rotation (about l/10th of 
the solar surface shear) and only slightly distorted by intrinsic variability. 

This result is at odds with various theoretical predictions, suggesting that dynamo 
fields of fully-convective stars should be mostly non- axisymmetric unless they succeed 
at triggering significant differential rotation. 

Key words: stars: magnetic fields - stars: low-mass - stars: rotation - stars: indi- 
vidual: V374 Peg - techniques: spectropolarimetry 
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1 INTRODUCTION 

Most cool stars exhibit demonstrations of activity such as 
cool surface spots, coronal activity or flares, usually at- 
tributed to magnetic fields. The current understanding is 
that these magnetic fields are generated at the surfaces 
and in the convective envelopes of cool stars through dy- 
namo processes, involving cyclonic motions of plasma and 
rotational shearing of internal layers. In partly-convective 
Sun- like stars, dynamo processes are thought to concentrate 
where rotation gradients are supposedly largest, i.e., at the 
interface layer between their radiative cores and convective 
envelopes. Cool sta rs with masses lower than 0.35 M be- 
ing fully convective ([Chabrier &; Bara ffe 1997), they lack the 
thin interface layer presumably hosting solar-type dynamo 
processes and are thought to rotate mainly as rigid bod- 
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ies (|Barnes et all 120051 : iKiiker Rudigerl l2005h : yet they 
are both very active (e.g.. jDelfosse et al.lll998|) and strongly 
magnetic (|Saar fc Linskvlll985l : IJohns-Krull fc Valentil[l996l : 
iReiners fc Basrill2007h 

This led theoreticians to propose that the intense ac- 
tivity and strong magnetism of fully convective dwarfs may 
be due to non-solar dynamo processes, in which cyclonic 
convection and turbulence play the main ro les while dif- 
feren tial rotation contributes very little (e.g., [Purnev et all 
1993); the exact mechanism capable of producing such fields 
is however still a debated issue. The most recent numerical 
dynamo simulations suggest that these stars are actually ca- 
pable of producing large-scale magnetic fields, although they 
disagree on the actual properties of the large scale field. 
While some find that fully- convective stars rotate as solid 
bodies and host purely non-axisymmetric large- scale fields 
(jKiiker fe Rudigerll2005l : IChabrier fe Kukerll2QQ6l ) , others di- 
agnose that these stars should succeed at triggerring differ- 
ential rota tion and thus produ ce a net axisymmetric poloidal 
field fe.g.. lDobler et alJl2006h . 

Existing observational data on fully convective dwarfs 
do not completely agree with any of these models. Observa- 
tions of low-mass stars indicate that surface differential rota- 
tion vanishes with increasing convective depth, so that fully- 
convective stars rotate mostly as solid bodies ([Barnes et al.l 
2005). At first sight, this result confirms nicely the predic- 
tions of lKiiker &; Riidigerl ((2005) , leading one to expect fully- 
conve ctive stars to host purely non-axisymmetric large-scale 
fields (jChabrier fe Kuke J2OO6I ) . However, spectropolarimet- 
ric data of the rapidly rotating M4 dwarf V374 Peg (G 188- 
38, GJ 4247, HIP 108706) collected with ESPaDOnS for al- 
most tln^_£om2^te_£otation periods over nine rotation cy- 
cles (jDonati et al.ll20Q6l hereafter D06) demonstrated that 
this fully- convective star hosts a strong mostly-axisymmetric 
poloidal field despite rotating almost rigidly, in contradiction 
with theoretical expectations. 

In this new paper, we present and analyse spectropo- 
larimetric data of V374 Peg collected with ESPaDOnS at 
three different epochs, 2005 August (i.e., those used in D06), 
2005 September (providing only partial coverage of the ro- 
tation cycle) and 2006 August (providing full, though not 
very dense, coverage of the rotation cycle). First, we briefly 
describe the new observations and the Zeeman detections 
we obtained. We then apply tomographic imaging to our 
data sets and characterise the spot distributions and mag- 
netic topologies at the surface of V374 Peg in both 2005 and 
2006; from the observed temporal evolution of the magnetic 
topology, we also estimate both the differential rotation and 
the long-term magnetic stability of V374 Peg. Finally, we 
discuss the implications of our results for our understanding 
of the dynamo processes operating in fully- convective stars. 



2 OBSERVATIONS 

Spectropolarimetric observations of V374 Peg were collected 
in 2005 August (hereafter Aug05), 2005 September (here- 
after Sep05) and 2006 August ( hereafter AugQ6) w it h the 
ESPaDOnS spectropolarimeter (|Donati et al.l Il997l I2006L 
Donati et al, in preparation) at the Canada-France-Hawaii 
Telescope (CFHT). The ESPaDOnS spectra span the entire 
optical domain (from 370 to 1000 nm) at a resolving power of 



about 65,000. During the Aug05 and Aug06 runs, V374 Peg 
was observed several times per night to ensure correct sam- 
pling of the rotational cycle (0.4456 d, D06). A total of 91 
circular-polarization (Stokes V) observations were collected, 
each consisting of four individual subexposures taken in dif- 
ferent polarimeter configurations (to suppress all spu rious 
polarisation signatures at first order, [Donati et aDll997h . 

Data reduction was performed using Libre ESpRIT 
(|Donati et al.l fl997l Donati et al 2007, in preparation), a 
fully automatic reduction package installed at CFHT and 
achieving optimal extraction of ESPaDOnS spectra. The 
peak signal-to-noise ratio (S/N) per 2.6 kms -1 velocity bin 
in the recorded spectra ranges from 99 to 222, depending 
mostly on weather conditions and exposure time (see Ta- 
ble [TJ. The Aug05 and Aug06 data sets respectively include 
64 and 21 spectra; 4 spectra in the Aug05 data set are signif- 
icantly affected by flaring and were discarded for the present 
study. Along with the 6 spectra collected in Sep05, the com- 
plete data set used for this study includes 87 spectra. 

In order to increase further the quality of our data set, 
we u sed Least-Squares Deconvolution (LSD, iDonati et al.l 
1 19971 ) to extract the polarimetric information from most 
photospheric spectral lines. The line list required for LSD 
was computed with an Atlas9 model atmosphere (|Kurucj 
1993) matching the properties of V374 Peg and includes 
about 5 000 spectral features. We finally obtain LSD Stokes 
V profiles with relative noise levels of about 0.05% (in units 
of the unpolarised continuum), corresponding to a typical 
multiplex gain of about ICQ with respect to individual spec- 
tra. A typical example is shown in Fig. [1] 

Zeeman signatures are detected in nearly all the spectra, 
with amplitudes of up to 0.5% (see Fig. [TJ and exhibit clear 
temporal variations. The red ward migration of the detected 
Zeeman signatures is particularly obvious in the Aug05 data 
set (see Sec. l4.2[) , demonstrating unambiguously that the ob- 
served variability is due to rotational modulation. Stokes / 
profiles only show a low (though clear) level of variability, 
with photospheric cool spots causing moderate line distor- 
sions (see Sec. 13. ip . 

In the following, all data are phased according to the 
following ephemeris (shifted by 0.01 phase with respect to 
that of D06): 

HJD = 2453601.78613 + 0.4456£. (1) 



3 MODELLING INTENSITY PROFILES 
3.1 Model description 

We use Doppler Imaging (DI) to convert time-resolved series 
of Stokes I LSD profiles into maps of cool spots at the sur- 
face of V374 Peg. Thanks to the Doppler effect, cool spots 
on the photosphere of a rapidly rotating star produce profile 



1 This multiplex gain refers to the peak S/N achieved in the 
spectrum, i.e., around 850 nm in the case of M dwarfs and ES- 
PaDOnS spectra. Previous studies, reporting higher gains (e. g. 
IDonati et al .11 19971 ), were based on spectra at shorter wavelengths 
(ranging from about 500 to 700 nm), where the line density at 
peak S/N is comparatively higher than in the present case. 
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Table 1. Journal of observations. Subsets of observations collected sequentially on the same night (with no time gaps) are listed on a 
single line. Columns 1-7 list the UT date, the number of observations in each subset, the heliocentric Julian date, the UT time, the 
exposure time, the peak signal to noise ratio (per 2.6 kms -1 velocity bin) and the rms noise level (relative to the unpolarised continuum 
level and per 1.8 kms -1 velocity bin) in the average circular polarization profile produced by Least-Squares Deconvolution (see text). 
The rotational cycle E from the ephemeris of equation (pQ) is given in column 8. Whenever several observations are present in a subset, 
cols. 3, 4 and 8 list the values corresponding to the first and last subset observations, whereas cols. 6 and 7 mention the minimum and 
maximum values reached within the subset. 



Date 


^obs 


HJD 


UT 




t 


6Xp 


S/N 


CLSD 


Cycle 






(2,453,000+) 


(h:m:s) 




(s) 




(io- 


4 /c) 




2005 


Aug 19 


6 


601.78613-601.87607 


06:43:29-08:52:59 


4 


X 


300.0 


145-170 


6.0- 


-7.7 


0.00-0.20 


Aug 19 


12 


601.94413-602.13790 


10:30:60-15:10:01 


4 


X 


300.0 


99-172 


5.6- 


11.3 


0.35-0.79 


Aug 21 


23 


603.75042-604.19199 


05:52:00-15:15:52 


4 


X 


300.0 


137-208 


5.7- 


-7.8 


4.41-5.40 


Aug 23 


23 


605.73892-606.13069 


05:35:25-14:59:33 


4 


X 


300.0 


125-176 


6.0- 


-8.8 


8.87-9.75 


Sep 08 


6 


621.94402-622.02801 


10:30:49-12:31:46 


4 


X 


300.0 


135-161 


7.1- 


-9.2 


45.24-45.43 


2006 


Aug 05 


2 


952.92084-952.93771 


09:58:04-10:22:22 


4 


X 


275.0 


172-181 


5.7- 


-5.8 


788.00-788.04 


Aug 05 


1 


953.05250 


13:07:39 


4 


X 


275.0 


180 


5.6 


788.30 


Aug 07 


2 


955.02780-955.04651 


12:31:58-12:58:55 


4 


X 


300.0 


192-192 


5.3- 


-5.3 


792.73-792.78 


Aug 08 


2 


955.86045-955.87780 


08:30:57-08:55:55 


4 


X 


300.0 


163-182 


5.5- 


-6.4 


794.60-794.64 


Aug 08 


2 


956.02162-956.03903 


12:23:01-12:48:05 


4 


X 


300.0 


194-198 


5.2- 


-5.4 


794.96-795.00 


Aug 09 


2 


956.82353-956.84159 


07:37:44-08:03:44 


4 


X 


300.0 


185-191 


5.2- 


-5.5 


796.76-796.80 


Aug 09 


2 


957.02738-957.04305 


12:31:15-12:53:50 


4 


X 


260.0 


182-187 


5.3- 


-5.4 


797.22-797.26 


Aug 10 


2 


957.84510-957.86535 


08:08:44-08:37:54 


4 


X 


360.0 


212-222 


4.6- 


-4.6 


799.06-799.10 


Aug 10 


2 


958.03059-958.05071 


12:35:50-13:04:49 


4 


X 


360.0 


191-208 


4.7- 


-5.1 


799.47-799.52 


Aug 11 


2 


959.03187-959.05226 


12:37:38-13:06:59 


4 


X 


360.0 


215-220 


4.7- 


-4.7 


801.72-801.76 


Aug 12 


2 


960.03589-960.05609 


12:43:23-13:12:27 


4 


X 


360.0 


205-205 


4.8- 


-4.9 


803.97-804.02 



V374 Peg, LSD profiles 
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Figure 1. Typical examples of LSD Stokes V and I profiles of 
V374 Peg (top and bottom curves respectively) on 2006 August 08 
(green line) and 2006 August 10 (black line). 



distorsions in spectral lines whose locations strongly corre- 
late with the spatial positions of the parent spots; in this 
respect, spectral lines of a spotted star can be viewed as one- 
dimensional (ID) image, resolved in the direction perpendic- 
ular to both the stellar rotation axis and the line of sight. By 
looking at how these ID maps are periodically modulated as 
the star rotates enables one to recover a 2D map of the stel- 
lar photosphere fe.g., IVogt et al.lll987h . Since this inversion 
problem is partly ill-posed, one needs to implement addi- 



tional constraints to stabilise the imaging process, e.g., by 
selecting the image having the lowest inf ormation content 
using the maximum entropy algorithm of ISkilling fe Brvanl 
(1984). More detai ls about the princ iples and performance 
DI can be found in lVogt et all (| 19871 ). 

We implement this imaging process by dividing the sur- 
face of the star into a grid of 3,000 elementary cells. We 
describe the local line profile at each grid point of the ste l- 
lar surface with the 2-component model of lCamero n (1992), 
i.e., as a linear combination between two reference profiles, 
one representing the quiet photosphere (at temperature T q ) 
and one representing cool spots (at temperature T s ), both 
Doppler shifted by the local line-of-sight rotation velocity 
of the corresponding grid cell. The image quantity we re- 
construct is the local spottedness at the surface of the star, 
i.e., the amount by which both reference profiles are mixed 
together, varying between (no spot) and 1 (spots only). 
The algorithm aims at reconstructing the image with mini- 
mum spot coverage at the surface of the star, given a certain 
quality of the fit to the data. 

More specifically, we assume that both reference profiles 
are equal and only differ by their relative continuum levels 
(set to their respective blackbody fluxes at the mean line 
wavelength). We further assume T q = 3200 K (i.e., a typical 
surface temperature for a M4 star) and T s = 2800 K (i.e., 
a low spot-to-photosphere temperature cont rast in agree- 
ment with the findings of iBerdvuginal 12005). For the as- 
sumed reference profile, we can use either the LSD pro- 
file of the very sl owly rotating inactive M dwarf Gl 402 
(jDelfosse et alJ ll998), or a simple gaussian profile with simi- 
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lar full-width-at-half-maximum (FWHM, set to 9.25 kms -1 
including instrumental broadening) and equivalent width. 
Both options yield very similar results, demonstrating that 
the exact shape of the assumed local profile has very minor 
impact on the reconstructed images, as long as the projected 
rotation velocity of the star is much larger than the local 
profile. We further assume that continuum limb- darkening 
varies linearly with the cosi ne of the limb angle (with a slope 
of u = 0.6964, IClaretll2004r ); using a quadratic (rather than 
linear) dependence produces no visible change in the result. 

The angle i of the stellar rotation axis to the line-of- 
sight is taken to be 70° (e. g., D06). Given the projected 
equatorial velocity vsini of V374 Peg (about 36.5 kms -1 , 
see below) and the width of the local profile (9.25 kms -1 , 
see above) , the number of spatially resolved elements across 
the equator of the star is about 25, which is equivalent to 
a longitude resolution of 15° or 0.04 rotation cycle. Using 
3,000 grid cells at the surface of the star (112 across the 
equator) is therefore perfectly adequate for our needs. 



line profiles are rather small. Note that only spots (or spot 
groups) with sizes comparable to or larger than the reso- 
lution element can be reconstructed with DI, tiny isolated 
spots much smaller than our spatial resolution threshold 
remaining out of reach. The resulting brightness maps are 
shown in Fig. 2] While both maps are roughly similar, they 
apparently differ by more than a simple phase shift (that 
could result from a small error in the rotation period). Note 
that no obvious polar spot is reconstructed on V374 Peg. 
This result is quite reliable; since V374 Peg is a fast rotator, 
any strong polar structure would produce a clear distortion 
in the core of the Stokes I profiles, giving them an obvious 
flat-bottom shape. This is not the case here, as one can see 
on Fig. [2] 

Small changes in vsini (within 0.5 kms -1 ), in the in- 
clination angle (within 10°), in the local profile model (see 
above) or in the limb- darkening model (linear vs quadratic) 
produce negligible modifications in the reconstructed spot 
distributions. 



3.2 Results 

DI is very sensitive to the assumed vsini (e.g., IVogt et aD 
Il987h : small errors are known to generate specific belt-like 
artifacts around the stellar equator in the reconstructed im- 
age. We derive the optimum v sini by minimising the recon- 
structed spot coverage at a given reduced chi-square (here- 
after Xr) ^ to the data, and find values of 36.5 ±0.4 kms -1 
and 36.7 =b 0.4 kms -1 for the Aug05 and Aug06 data sets 
respectively, when using the gaussian local profile model. 
Slightly different values (all lying with 0.5 kms -1 of the pre- 
vious estimates) are obtained when using the LSD profile of 
Gl 402 as local profile, or when using different models of the 
continuum limb darkening. Taking into account all potential 
sources of systematic errors, we obtain that the absolute ac- 
curacy to which vsini is determined is about 1 kms -1 ; in 
the following, we assume vsini = 36.5 kms -1 . 

Aug05 and Aug06 datasets can be fitted at S/N lev- 
els of 715 and 800 respectively. In particular, the distorsions 
travelling from the blue to the red profile wing are well repro- 
duced (see Fig. [2] and [3}. We note small systematic residual 
discrepancies between the observed and modelled profiles, 
showing up as vertical bands in the far line wings (see Fig. [2] 
and [3]). These differences cannot be removed by adjusting 
vsini or u and do not affect the imaging of starspots. An- 
other residual discrepancy is visible in the blue profile wing 
around phase 0.27 in the Aug06 data set (see Fig. [3] right 
panel). This discrepancy corresponds to a feature visible in 
the blue wing of the observed profiles (see Fig. [3] left panel), 
that the imaging code did apparently not convert into a 
starspot. The reason for this is that the code assumes the 
rotating star to host starspots exhibiting no intrinsic vari- 
ability; if the observed (unfitted) feature was due to such 
a starspot, we would expect a similar signature to show up 
in the red wing of the observed profiles around phase 0.5, 
i.e., when the putative spot has reached the receding stellar 
limb. Since this is not the case here, we speculate that the 
residual unfitted spectral feature is either caused by a small 
flare or by an intrinsically variable starspot. 

Both datasets yield spot occupancies of about 2% of the 
overall photosphere at both epochs. This rather low spotted- 
ness directly reflects the fact that distortions of unpolarised 



4 MODELLING CIRCULARLY POLARISED 
PROFILES 

4.1 Model description 

We use Zeeman-Doppler Imaging (ZDI) to turn our series 
of rotationally modulated circularly-polarised LSD profiles 
into magnetic field maps at the surface of V374 Peg. ZDI is 
based on the same basic principles as DI except that it aims 
at mapping a vector field. The overall characteristics of the 
model we use are given below. 

The vector magnetic field at the surface of the star is de- 
scribed using spherical harmonics expansions f or each of the 
field c ompo nent in spherical coordinates (e.g., [Donati et aD 
l200ll . l2006h . Note that this is different from the initia l 
ZDI attempts (e.g., iBrown et aDll99ll : iDonati et aDll997h 
in which the field was described as a set of independent 
pixels, containing the components of the field in each sur- 
face grid cell. The main advantage of this new method is 
that it straightforwardly allows the recovery of a physically 
meaningfull magnetic field, that we can easily split into its 
poloidal and toroidal components. This is of obvious inter- 
est for all studies on stellar dynamos. Another important 
point is that this new method is very successfull at recov- 
ering simple magnetic field structures such as dipoles (even 
from Stokes V data sets only), while the old one reportedl y 
failed at such tasks ([Brown et al.lll99ll : iDonati et al.ll200lh . 

As for DI (see Sec. 13. 1[) , the stellar surface is divided 
into a grid of 3 000 pixels on which the magnetic field com- 
ponents are computed from the coefficients of the spherical 
harmonics expansion. The contribution of each grid point to 
the synthetic Stokes V spectra are derived using the weak- 
field approximation: 

where / is the local unpolarised profile (modelled with the 
same gaussian profile as in Sec. 13. ip . g and Ao the Lande 
factor and the centre rest wavelength of the average LSD 
line (set to 1.2 and 700 nm respectively), e and m the elec- 
tron charge and mass respectively, and Bi the line-of-sight 
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Figure 2. Observed (left), modelled (middle) and residual (right) Stokes I dynamic spectra of V374 Peg for the Aug05 data set (60 
spectra). A purely rotational profile is subtracted from all spectra to emphasize the image contrast. The fit corresponds to a S/N level 
of 700. The colour scale ranges from -0.5% to 0.5% of the continuum level, red and blue respectively standing for a lack and an excess 
of absorption with respect to the synthetic purely rotationnal profile. 



All 



II 



I 



n 



■I I 



lh 



UHU.I 



J 



■JJ 



i ■ hi i 

h ■ ■■ ■ k 



-50 

Velocity (km/s) 




J_i 



M 



1 1 



nvT i 

ii i in i 

I ■!! I 



¥ v 



II 1 1 II II I 

ft i>U 



Velocity (km/s) 



-50 50 

Velocity (km/s) 



Figure 3. Same as Fig. [2]for the Aug06 data set (21 spectra). 



projection of the local magnetic field vector. Despite being 
an approximation, this expression is valid for magnetic split- 
tings significantly smaller than the local profile width (prior 
to instrumental broadening); given that 2 kG fields yield 
Zeeman splittings of only about 2.5 kms -1 in the case of 
our average LSD profile (i.e., about twice smaller than the 
FWHM of the local profile prior to instrumental broaden- 
ing), this expression turns out to be accurate enough for 
our needs. This is further confirmed by using model local 



profiles computed from Unno-Rachkovsky's equations, with 
which almost identical results are obtained. 

As for DI, the inversion problem is partly ill-posed. We 
make the solution unique by introducing an entropy function 
describing the amount of information in the reconstructed 
image, and by looking for the image with maximum entropy 
(minimum information) among all those fitting the data at a 
given Xr level. The entropy is now computed from the sets of 
complex spherical harmonics coefficients, rather than from 
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Figure 4. Doppler maps of the spottedness at the surface of V374 Peg reconstruct, as derived from the Aug05 (left panel) and Aug06 
(middle panel) data sets. The last image (right) was recovered from the complete dataset (Aug05, Sep05, and Aug06 altogether), it 
represents the starspots configuration at an intermediate epoch (cycle 194) and takes into account the adjustment of the period and the 
inclusion of differential rotation (see Sec. 5). The star is shown in flattened polar projection down to latitudes of —30°, with the equator 
depicted as a bold circle and parallels as dashed circles. Radial ticks around each plot indicate phases of observations. 



the individual image pixels as before. We chose one of the 
simplest possible form for the entropy: 

S = - (a|, m + film + llm) (3) 

where a^, m , /3g,m, j£,m are the spherical harmonics coeffi- 
cient of order (£, m) describing respectively the radial, non- 
radial poloidal and toroidal field c omponents. More d etail 
about the process can be found in iDonati et ail (|2006T ). In 
particular, this function allows for negative values of the 
magnetic field (as opposed to the conventional expression of 
the Shannon entropy). The precise form of entropy (e.g., the 
multiplying factor £) has only minor impact on the resulting 
image. 

As mentioned in in Sec. 13.11 the number of resolved el- 
ements on the equator is about 25, indicating that limiting 
the spherical harmonics expansion describing the magnetic 
field to orders lower than 25 should be enough and gener- 
ate no loss of information. In practice, no significant change 
of the reconstructed image is observed when limiting the 
expansion at orders < 10. All results presented below are 
obtained assuming £ < 10. 

4.2 Results 

We present here the results concerning the Aug05 dataset, 
though already reported in D06, in order to clarify the com- 
parison with the Aug06 observations. The Aug05 data set 
can be fitted with a purely poloidal field down to %r = 1-1 
(see Fig. [5] for the observed and modelled dynamic spectra 
at this epoch). No improvement to the fit is obtained when 
assuming that the field also features a toroidal component. 
We therefore conclude that the field of V374 Peg is mostly 
potential in Aug05 and that the toroidal component (if any) 
includes less than 4% of the overall large-scale magnetic field 
energy. The field we recover is shown in Fig. (left column) 
while the corresponding spherical harmonics power spectra 
are shown in Fig. [8] (left column). The reconstructed large- 
scale field is found to be mostly axisymmetric, with the vis- 
ible hemisphere being mostly covered with positive radial 
field (i.e., with radial field lines emerging from the surface 
of the star). As obvious from Fig. [8] the dominant mode 
we recover corresponds to £ = 1 and m — (i.e., a dipole 
field aligned with the rotation axis) and has an amplitude 



of about 2 kG; all other modes have amplitudes lower than 
1 kG. The maximum field strength at the surface of the star 
reaches 1.3 kG, while the average field strength is about 
0.8 kG. 

Similar results are obtained from the Aug06 data, for 
which the observed spectra are fitted down to %^ = 1.0 (see 
Fig. [6}. Again, we find that the field is mostly axisymmetric, 
with 80% of the reconstructed magnetic energy concentrat- 
ing in m = modes. The radial field component clearly 
dominates (75% of the energy content) while the toroidal 
component is once more very marginal (4% of the energy 
content). The main excited mode again corresponds to £ = 1 
and m — 0, features an amplitude of 1.9 kG and encloses 
as much as 60% of the overall magnetic energy content. The 
maximum and average field strengths recovered at the sur- 
face of the star are 1.2 kG and 0.7 kG respectively. 

As a result of the limited spatial resolution of the imag- 
ing process, small-scale magnetic regions at the surface of 
V374 Peg are out of reach of ZDI. We can therefore not eval- 
uate the amount of magnetic energy stored at small scales 
and compare it to that we recover at medium and large 
scales. We can however claim that scales corresponding to 
spherical harmonics orders ranging from about 5 to 25, to 
which our imaging process is sensitive, contain very little 
magnetic energy with respect to the largest scales (£ < 5). 

We can also note that that the magnetic maps recov- 
ered from the Aug05 and Aug06 data sets are fairly similar 
apart from a 36° (0.1 rotation cycle) global shift between 
both epochs, certainly resulting from a small error in the 
rotation period. Slight intrinsic differences also seem to be 
present, e.g., in the exact shape of the positive radial field 
distribution over the visible hemisphere. A straightforward 
visual comparison of both images is however not sufficient 
to distinguish whether these apparent discrepancies simply 
result from differences in phase coverage at both epochs, 
or truly reveal intrinsic variability of the magnetic topol- 
ogy between our two observing runs. To test this, one way 
is to merge all available data (including those collected in 
Sep05) in a single data set and to attempt fitting them with 
a single magnetic structure (and a revised rotation period) . 
This is done in Sec. [5] In the meantime, we can nevertheless 
conclude that the magnetic topology of V374 Peg remained 
globally stable over a timescale of 1 yr. 
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Figure 5. Observed (left), modelled (middle) and residual (right) Stokes V dynamic spectra of V374 Peg for the Aug05 data set (60 
spectra). The fit corresponds to Xr = 1-1- The colour scale ranges from -0.4% to 0.4% of the continuum level, red and blue respectively 
standing for positive and negative circular polarization signal. 
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Figure 6. Same as Fig. for the Aug06 data set (21 spectra). The fit corresponds to Xr — 1-0- 



5 DIFFERENTIAL ROTATION AND 

MAGNETIC FIELD TIME-STABILITY 

To estimate the degree at which the magnetic topology re- 
mained stable over 1 yr, we merge all data together (Aug05, 
Sep05 and Aug06 datasets) and try to fit them simultane- 
ously with a single field structure. We also assume that the 
surface of the star is rotating differentially (with a given dif- 
ferential rotation law, see below) over the complete period 
of our observations, and thus that the magnetic topology 



is slightly distorted as a result of this surface shear (with 
respect to the magnetic topology at median epoch, e. g. at 
cycle 194). If the fit to the complete data set is as good as 
that to both individual Aug05 and Aug06 data sets (xr — lj 
see Sec. I4.2[) , it would demonstrate that no variability oc- 
curred (except that induced by differential rotation) over 
the full observing period. If, on the other hand, the fit qual- 
ity to the complete data set is worse, it would indicate that 
intrinsic variability truly occurred, with the fit degradation 
directly informing on the actual level of variability. 
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Figure 7. Surface magnetic field of V374 Peg as derived from our Aug05 (left) and Aug06 (middle) datasets. The last image (right) was 
recovered from the complete dataset (Aug05, Sep05, and Aug06 altogether), it represents the magnetic configuration at an intermediate 
epoch (cycle 194) and takes into account the adjustment of the period and the inclusion of differential rotation (see Sec. 5). The three 
components of the field in spherical coordinates are displayed from top to bottom (flux values labelled in G). The star is shown in 
flattened polar projection down to latitudes of —30°, with the equator depicted as a bold circle and parallels as dashed circles. Radial 
ticks around each plot indicate phases of observations. 



To model d iffere ntial rotation , we proceed as in 
IPetit et all (j2002h and iDonati et all (|2003h . assuming that 
the rotation rate at the surface of the star is solar- like, i.e., 
that Q varies with latitude as: 

Q(0) = Q eq -dQsm 2 (4) 

where Q eq is the rotation rate at the equator and the 
difference in rotation rate between the equator and the pole. 
We use this law to work out the phase shift of each surface 
grid cell, at any given epoch, with respect to its position at 
median epoch. These phase shifts are taken into account to 
compute the synthetic spectra corresponding to the current 
magnetic field distribution at all observing epochs using the 
model described in Sec. 14.11 

For each pair of Q eq and dQ values within a range 



of acceptable values (given by D06), we then derive, from 
the complete data set, the corresponding magnetic topol- 
ogy (at a given information content) and the associated 
Xr level at which modelled spectra fit observations. By fit- 
ting a paraboloid t o the Xr surface derived in this process 
((Donati et al.ll2003l ). we can easily infer the magnetic topol- 
ogy that yields the best fit to the data along with the cor- 
responding differential rotation parameters and error bars. 

The best fit we obtain corresponds to Xr — 1-15, demon- 
strating that intrinsic variability between Aug05 and Aug06 
remains very limited, as guessed in Sec. 14.21 by visually 
comparing the images from individual epochs. The derived 
Xr m &p (see Fig. [5} features a clear paraboloid around 
the x 2 minimum, yielding differential rotation parameters 
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Figure 8. Spherical-harmonics power spectra of the magnetic 
field maps derived from our Aug05 (left), Aug06 (middle) and 
complete (right) data sets (magnetic fluxes labelled in G) after 
adjustment of the period and inclusion of differential rotation (see 
Sec.0. From top to bottom, the magnetic power spectra respec- 
tively correspond to the radial poloidal, the non-radial poloidal 
and the toroidal field components. Only modes up to order 5 are 
plotted here (larger order modes containing very little power). 




1.153 



14.0984 



14.0992 



■1.142 



Figure 9. Variation of Xr as a function of the differential rotation 
parameters f2 e q and dfl. The last colour contour corresponds to 
a Xr increase of 1%, i.e., to an error bar of about 7a (for each 
parameter taken separately). 



range of latitudes (as r equested for a cle an differential rota- 
tion measurement, see Petit et al.l [20021 ) but rather concen- 
trate at low latitudes. Using the differential rotation param- 
eters derived from Stokes V data yields a fit corresponding 
to a S/N level of about 680, i.e., only slightly lower than 
that achieved at each individual epoch (see Sec. I3.2[) . It sug- 
gests that only moderate intrinsic variability occurred in the 
brightness distribution as well, despite the apparent (but 
weakly significant) differences between the two individual 
images. The corresponding image (reconstructed from the 
full data set) is shown in Fig. U (right panel). 



at the surface of V374 Peg equal to Q eq = 14.09880 ± 
0.00006 radd -1 (corresponding to P ro t = 0.445654 ± 
0.000002 d) and = 0.0063 ± 0.0004 radd -1 . Due to 
the large time gap between the two main observing epochs 
(about lyr), the Xr map exhibits aliases on both left and 
right sides of the minimum, corresponding to shifts of ^ 
6 x 10 -4 d on Prot- The nearest local minima, located at 
Q eq = 14.08110 radd -1 and Q eq = 14.11664 radd" 1 , are as- 
sociated with Xr values of 1.27 and 1.30 respectively, i.e., to 
Ax 2 values of 600 and 740 respectively; the corresponding 
rotation rates can thus be safely excluded. 

In addition to refining the rotation period (error bar 
shrunk by 2 orders of magnitude with respect to D06), we 
also detect a small but definite surface shear, about 10 times 
weaker than that of the Sun. Assuming solid body rotation 
does not allow to reach Xr levels below 1.5; this unambigu- 
ously demonstrates that the surface of V374 Peg is not ro- 
tating as a solid body and is truly twisted by differential 
rotation. The reconstructed magnetic map along with the 
corresponding spherical-harmonics power spectra are shown 
in the right panels of Figs. [7] and [8] 

The same procedure is used on Stokes I data. However, 
the Xr m &P we obtain does not feature a clear paraboloid, 
but rather a long valley with no spatially well-defined min- 
imum; this is not unexpected given the fact that spots re- 
constructed at the surface of V374 Peg do not span a wide 



6 CHROMOSPHERIC ACTIVITY 

We compute the mean and standard deviation profiles of 
the Ha, K/3 and Ca n infrared triplet (IRT) lines. Both Ha 
and H/3 lines are in strong emission whereas the infrared 
triplet is observed mostly in absorption with sporadic emis- 
sion episodes (due to flares). We notice that, for all studied 
lines, variability concentrates within ivsini about the line 
centre (see Fig. 10), suggesting that activity is mostly lo- 
cated close to the surface. 

For each observation, we compute the emission equiva- 
lent width of the chromospheric lines. The result for Ha is 
plotted on Fig. 11 as a function of rotational phase, H/3 and 
Ca II IRT exhibiting a very similar behaviour. We first con- 
clude that all three lines do not undergo pure rotational 
modulation, sometimes displaying very different emission 
levels at similar phases of successive rotation cycles (e.g., 
phase 0.2, see Fig. 11). We attribute the basal emission com- 
ponent (in Balmer lines) to quiescent activity and the intrin- 
sically variable emission component to flares. 

Maximum flaring activity is observed around phases 
0.20 and 0.75-1.00, but also, to a lesser extend at phases 
0.50-0.60 (see Fig. 11, bottom panel). Note that the four 
photospheric profiles discarded from our imaging analysis 
(see Sec. [2} were collected on 2005 August 21 at phases 0.72 
and 0.78 and on 2005 August 23 at phases 0.19 and 0.23, 
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Figure 10. Ha mean profile (black solid line) and standard de- 
viation (red solid line, shifted by +1), along with continuum level 
(dashed line). All profiles are shown in the stellar velocity rest 
frame. 




Rotationnal phase 



Figure 11. Emission equivalent width of Ha as a function of 
the rotational phase in Aug05 and Sep05, consecutive observa- 
tions are linked together by a solid line (top panel); along with 
associated standard deviation (computed over phase intervals of 
0.1 cycle, bottom panel). 



i.e., at the very beginning of the two main flaring episodes 
we recorded. 



7 DISCUSSION AND CONCLUSION 

Spectropolarimetric observations of V374 Peg were carried 
out with ESPaDOnS and CFHT at 3 epochs over a full time 
span of about 1 yr. Clear Zeeman signatures are detected in 
most Stokes V spectra; weak distortions are also observed 
in the unpolarised profile of photospheric lines. 

Using DI, we find that only a few low-contrast spots are 
present at the surface of V374 Peg, covering altogether 2% 
of the overall surface. The spot distribution is apparently 
variable on a timescale of 1 yr. Even though we may be 
missing small spot s evenly peppering the photosphere (e.g., 
Ijeffers et al.1 12006), V374 Peg nevertheless appears as fairly 



different from warmer, partly-convective, active stars, whose 
surfaces generally host large, high-contrast cool spots often 
covering most of the pola r regions at rotatio n rates as high 
as tha t of V374 Peg fe.g.. iDonati et alllioooh . iBarnes et all 
(2004) also report the lack of a cool polar structure for the 
Ml. 5 - largely but not fully convective - fast rotator HK Aqr 
(P ro t = 0.431 d). A spot occupancy of about 2% has been 
observed for very s low rotators as well (|Bonfils et al.ll2007l : 
Pe rnor v et al. I l2007h . suggesting that spot coverage may not 
correlate with rotation rate in M dw arfs, as op pposed to 
what is observed in hotter stars (e. g., lHall|[l99lh . 

Using ZDI, we find that the surface magnetic topology 
of V374 Peg is mostly poloidal and fairly simple, with most 
of the energy concentrating within the lowest-order axisym- 
metric modes. This confirms and amplifies the previous re- 
sults of D06. Again, this is at odds with surface magnetic 
topologies of warmer active stars, usually showing strong 
and often dominant toroidal field components in the form of 
azim uthal field rings mo re or less encircling the rotation axes 
(e.g. jDonati et al.ll2003h . As for surface spots, ZDI is mostly 
insensitive to small-scale magnetic regions (and in particular 
small-scale bipoles) that can also be present at the surface 
of V374 Peg; we can however claim that very little mag- 
netic energy concentrates in spherical-harmonic terms with 
£ ranging from 5 to 25. 

We find that the surface magnetic topologies of 
V374 Peg at epoch Aug05 and Aug06 are very similar apart 
from an overall shift of about 0.1 rotation cycle. This is con- 
firmed by our success at fitting the complete Stokes V data 
set down to almost noise level with a unique magnetic field 
topology twisted by a very small amount of solar-like dif- 
ferential rotation. It suggests in particular that magnetic 
topologies of fully- convective stars are globally stable on 
timescales of at least 1 yr. Again, this is highly unusual 
compared to the case of warmer active stars, whose detailed 
magnetic topologies evolv e beyond recognitio n in timescales 
of only a few weeks (e.g., IDonati et al.ll2003l ). 

The updated rotation period and photospheric shear 
we find for V374 Peg are respectively equal to 0.445654 ± 
0.000002 d and 0.0063 ± 0.0004 rad d" 1 respectively; the cor- 
responding time for the equator to lap the pole by one com- 
plete rotation cycle is equal to 2.7 yr, about 10 times longer 
than solar. We speculate that this unusually small differen- 
tial rotation and long lap time are in direct relationship with 
the long lifetime of the magnetic topology, suggesting that 
differential rotation is what mostly controls the lifetime of 
magnetic topologies. New data collected on a smaller time 
span (a few months) and featuring in particular no intrin- 
sic variability (other than that generated by differential ro- 
tation) are needed to provide an independent and definite 
confirmation of this result. 

Although the spot distributions at the surface of 
V374 Peg at epoch Aug05 and Aug06 apparently differ from 
more than a simple phase shift of 0.1 cycle, both data sets are 
compatible with a unique brightness map undergoing mod- 
erate intrinsic variability (assuming differential rotation is 
similar to that derived from Stokes V data). The low spot- 
tedness level and weak spot contrast implies that densely 
sampled data sets at 2 nearby epochs are needed to confirm 
this result, derive an independent determination of differen- 
tial rotation parameters from unpolarised spectra and check 
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that spots and magnetic fields trace the same atmospheric 
layers. 

The strong mostly-axisymmetric poloidal field of 
V374 Peg can hardly be reconciled with the very weak 
level of surface differential rotation in the context of the 
most recent dynamo models of fully- convective stars. While 
Dobler et al.l (|20Q6h predicts that rapidly-rotating fully- 
convective stars should be able to sustain strong poloidal 
axisymmetric fields with the help of significant ( thoug h an- 
tisolar) differential rotation, IChabrier &; Kiikerl |2006) find 
that these stars should only be able to trigger purely non- 
axisymmetric magnetic topologies if differential rotation is 
lacking (or remains below a threshold of about 1/3 of the so- 
lar value, Chabrier, private communication). Both findings 
are at odds with our observational results, suggesting that 
our theoretical understanding of dynamo action in fully- 
convective stars is still incomplete. 

We report the detection of strong chromospheric ac- 
tivity in V374 Peg, including quiescent emission in Balmer 
lines along with sporadic flaring events in Ha, H/3 and the 
Call IRT. This is noticeably d i fferen t from what was re- 
cently reported by Berger et al . (2007) for the M8.5 brown 
dwarf TVLM 513-46546, whose rotationally modulated Ha 
emission identically repeats over several successive rotation 
cycles. From the width and variability of emission profiles, 
we conclude that activity occurs mostly close to the stellar 
surface on V374 Peg, with flares apparently concentrating at 
specific longitudes. However, there is no clear spatial correla- 
tion between the photospheric magnetic field, the starspots 
and these chromospheric regions, e.g., strongest magnetic 
fields being observed when chromospheric activity is low. 
Simultaneous mult i wavelength observations of V374 Peg, in- 
cluding in particular optical spectropolarimetry, X-ray, and 
radio monitoring, could allow us to investigate the 3D mag- 
netosphere of fully convective stars (e.g., by extrapolating 
the photospheric field topology up to the corona), and to 
unravel how photospheric activity relates to chromospheric 
and coronal activity. 

Finally, we note that our v sin i and rotation period es- 
timates imply that R* sini = 0.32 ± 0.01 R©, and thus that 
R* c± 0.34 R© (assuming i c± 70°). Given the Mr magni- 
tude of V374 Peg (equal to 7.0 =b 0.3), we ca n safely derive 
that V374 Peg has a mass of 0.28 ± 0.05 M (jPelfosse et al.l 
2000). It implies that the radius of V374 Peg is slightly larger 
than what the oretical mode ls predict for low-mass stars of 
similar masses. iRibasI (|2006h also reports stellar radii larger 
than the theoretical predictions for very-low mass eclips- 
ing binaries (whose rotation period is very short), whereas 
the radii of slowly rotating M dwa rfs match the predicted 
values (e. g.. ISegransan et a l. 2003). Ac cordin g to the phe- 
nomenological model of Chabrier et al.l ([2007) the presence 
of a strong magnetic field may inhibit stellar convection and 
therefore lead to larger radii with respect to an inactive star 
of equal mass. 

Future observations are thus needed to explore in more 
details the magnetic topologies, brightness distributions and 
differential rotation of V374 Peg and other fully- convective 
stars. Finding out how these quantities vary with mass 
and rotation rate could in particular show us the way to 
more successfull simulations of dynamo processes in low- 
mass fully-convective stars. 
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